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Improvement of the hot QCD pressure by the minimal sensitivity 

criterion 

Masahiro Inui*), Akira Niegawa**) and Hideaki Ozaki***) 
Graduate School of Science, Osaka City University, Osaka 558-8585, Japan 

The principles of minimal sensitivity (PMS) criterion is applied to the perturbative free 
energy density, or pressure, of hot QCD, which include the ~ lng a and part of the ~ 
terms. Applications are made separately to the short- and long-distance parts of the pressure. 
Comparison with the lattice results, at low temperatures, shows that the resultant "optimal" 
approximants are substantially improved when compared to the MS results. In particular, 
for the realistic case of three quark flavors, the "optimal" approximants are comparable with 
the lattice results. 

§1. Introduction 

In statistical physics, the free energy density, or pressure, plays a central role. 
At high temperature, the perturbative expansion of the pressure of the quark-gluon 
plasma is known to fifth order 1 ) in g s (the QCD coupling constant) and partially to 
sixth order. 2 ) These results are reliable only when g s is sufficiently small. Because 
of the asymptotic freedom of QCD, this is achieved at high temperature T > 2 GeV 
which is order of magnitude larger than the critical temperature (T c ) of the QCD 
phase transition. For T c < T < 2 GeV, the accuracy of the approximation, within 
the MS scheme, does not improve by the inclusion of higher order terms (conver- 
gence issue). Moreover, the unphysical dependence on the renormalization scale /i 
is strong. With the traditional choice [i = 2ttT, the truncated perturbation series 
for the pressure, in the MS scheme, does not reconcile with the lattice results (0.3 
GeV < T < 0.5 GeV). To get rid of these dilemmas, various kinds of resummations 
are proposed. Among those are quasi-particle models, 3 ) ^-derivative approxima- 
tions, 4 ) screened perturbation theory, 5 )' 6 ) hard thermal loop perturbation theory, 7 ) 
optimizations based on the principles of minimal sensitivity (PMS) and the fastest 
apparent convergence (FAC) criterions, 6 ) and the Pade type resummation and its 
variants. 8 )' 9 ) 

In this paper, putting aside the convergence issue, we will obtain''') the "optimal" 
approximants for the pressure p through an application of the PMS criterion. 10 ) 
Throughout this paper, for the expression p (in massless QCD), we use the one 
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Blaizot, Iancu, and Rebhan were the first who, among other things, applied the PMS criterion 
to the perturbative pressure. 6 ' The content of this (PMS) part of Ref. 6) and our way of improving 
it will be mentioned as occasions arise later. 
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given*) in Ref. 9): 

P = Pe + PM + PG , 

where pe, Pm, and pc come, in respective order, from the hard (~ T), soft (~ g s T), 
and supersoft (~ g 2 T) energy regions. The scales g s T and g 2 T are the color-electric 
and color-magnetic screening scales, respectively. A crucial observation made in Ref. 
9) is that the "short-distance pressure" pe and the "long-distance pressure" pm+g 
(= pm+Pg) ar e separately physical quantities, and thus both, when calculated to all 
orders of perturbation theory, are independent of the renormalization scheme. Then 
we will apply the PMS criterion to pe and pm+g separately**). 

In §2, a brief review of the PMS criterion is given. Then, we present the pertur- 
bation expansion of pe and pm+g m the form given in Ref. 9), which is taken as the 
basis of our analysis. In §3, we apply the PMS criterion to pe and pm+g separately 
and find the optimal approximants. In §4, we compare our results with the lattice 
results and the MS results. §5 is devoted to discussion and conclusion. 

§2. Preliminaries 

2.1. Brief review of the Principles of Minimal Sensitivity (PMS) (Ref. 10)) 

The couplant a s = g 2 / (47r) in massless QCD runs according to the renormaliza- 
tion group equation, 

^ = (3(a s ) = -a 2 s J2^i- (2-1) 
Here [i is the renormalization scale and 

(3i = ^ (51- y n/J , (2-2) 

with nf is number of quark flavors. The renormalization schemes can be labelled by 
a countable infinity of parameters (/i, •••)• [A) and Pi are independent of the 

renormalization scheme.] In the MS scheme, 

If 5033 325 9 \ , N 

(fc)ns = ^3 ( 2857 " — "/ + ^7 n f) ■ ( 2 - 3 ) 

The dependence of a s on (3j (j = 2, 3, ...) is governed by 
r)a f' as r J+2 

S = / 7^TT2 dx (= 7i(«.)) , (2-4) 



*' How to construct it from the available calculations are elucidated in Ref. 9) (see also Ref. 

6)). 

**' Since little information is available on pg, we will apply the PMS criterion to pm+g and not 
separately to pu and pc- Incidentally, in Ref. 6), the PMS criterion is applied to the total pressure 
p — Pe +Pm+g- 
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2n 

= g(pr^ a s +- (j = 2,3,...). (2-5) 

It is to be noted that a change in \i results in a change in a s of 0(a^) while a 
change in (3i {i = 2,3,...) results in a change in a s of 0(a* +1 ). The running a s = 
a s (fi, /?2, /?3, •••) is obtained by solving Eqs. (|2- 1|) and (|2-4j) under the boundary 
condition that is provided by the experiment. 

Let 7Z be some physical quantity, 1Z = a l s J2'jLo r j (X s- Then, when calculated to 
all orders of perturbation theory, 7Z is independent of the renormalization scheme, 
or, equivalently, of the parameters fx, 02, 03, ■■■ ■ 

The i-th order perturbative approximation 1Z® of 1Z is defined by (i) truncating 
the perturbation series for 1Z at ith order, and (ii) replacing the expansion parameter 
a s by its ith-order approximant defined as the solution to Eqs. (|2-1|) and (|2-4|) with 
P on the right-hand sides (RHS's) truncated at ith order; 

i-1 i-1 

K® =a l s ^2 r M > (<*s) = E Pi°i > ( 2 - 6 ) 

3=0 j=0 

= l? = 7fK), (2-7) 



where jj (a 8 ) 1S defined by Eq. (|2-4|) with 0® for 0. 

The PMS criterion 10 ^ for obtaining the optimal approximant for 1Z goes as fol- 
lows. 

1) Compute 0® and 1Z® in a particular renormalization scheme, say, MS scheme. 

2) Impose the requirement, 

In the following we refer Eqs. (|2-8|) and (|2-7|) to as the OPT equations. The 
OPT equations determine first (a) the dependence of the coefficients (r\, ...,rj_i) 
in Eq. ()2-6|) on (/•*> 02, ■ 0i-i), and then (b) the values for the parameters 
(/J-,02, — , A-l) = (^,02 > — >^i-l) together with a s = a£ . From (a) and (b), 
the values for the coefficients are determined, (ri, ...,rj_i) = (r^ , ...,r J -^ 1 ). 

3) The optimal approximant for TZ is obtained by substituting thus determined ats 

and (r^ , for, in respective order, a s and (n, ...,rj_i) in Eq. I|2-6|) . 

It should be noted here that the above definition of TZ® is not unique. Instead 
of (i) of the definition, one can adopt 

(i — l oo \ 

£ri«£ + 5>;°i ■ (2-9) 
i=o j=i / 

Here the coefficients (rj, rj_|_i, ...) are arbitrarily chosen fixed ones within the "rea- 
sonable" range. The original PMS is the special case, rj = rj+i = ... = 0, of this 
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more general setup. The PMS criterion can also be applied to more general case, 
e.g., 

oo rij 

K = J2 a * +j E ( r Jk + V^ s r' jk ) ln fc a s , (2-10) 

j=0 k=0 

where rij is or some integer. 

2.2. Quark- gluon plasma pressure 9 ) 

As has been mentioned in §1, the form for the pressure p given in Ref. 9) consists 
of three pieces, p = Pe+Pm+Pg- We will apply the PMS criterion to pe and pm+g 
(= Pm + Pg) separately. To define the border between the soft and hard energy 
regions, the scale parameter Ae is introduced, 0(g s T) < O(Ae) < 0(T). 

The long-distance pressure pm+g is given 9 ) as an expansion in powers of the 
effective EQCD parameters gE and the'- 



PM+G = —Tm%PM+G ■ 



Pm+g 



A e \ 110.2 



47T m E V 4 2m E J (47r) 2 



AX 

m E J 



81 



(A) L.441n. 
(4vr) 3 \m E J 

5 to, \ °% r 
- — (9 - nt) T, 

47r v m E 



A, 



2rriE 



+ 1.566 In 



A, 



1.392 + 5 G 



where 



m\ = 4vr(l + n f /6)a s T 2 



9% 



1 + a s [ (3 In 



1 + a 



ZTTl 7T 



1 



2,T + ^<"'> 



with 



Pm(nf) 



0.6124 - 0.4881n/ - 0.04280n 2 f 



(2-11) 



(2-12) 

(2-13) 
(2-14) 



l + n//6 

-0.3876 - 0.423571/. (2-15) 

The dimensionless number <5g in Eq. (j2-12|) is unknown. We will allow, as in Ref. 
9), the following variation: 

-5<(5 G <+5. (2-16) 
We refer 6 ) the expression (|2-11[) - (|2-12|) to as the untruncated form for pu+G- 

Expanding p^+Gi Eqs. (|2-lip and (|2-12p . in powers of a s and keeping up to 
and including O(af) contribution, we obtain the truncated form for pm+g, 



(2) _ (2) , (2) 

Pm+G ~ Pi, M+G ~r Pi, M+G 

Pi^m+g = Q s /2 ( r io + r na s ) ; 

P2?M+G = a l ( r 20 + r 2 l a 8 ) . 



(2-17) 
(2-18) 
(2-19) 
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We see from these expressions that, for pm+g> A 4 is the only renormalization-scheme 
parameter and p\ is not. The reason is as follows. On the one hand the highest-order 
term in p M \ G is of 0(p? s ) (see Eq. (|2-19)l ). On the other hand, a change in 02 results 

in a change in py, G of 0{aj 2 ) (< 0{a z s )) (cf. Eq. (|2-18|) and the comment above 
after Eq. (|2-5[Q . Thus, 02 is not the renormalization-scheme parameter. When the 
0(a s ) term of pm+g is computed, p\ becomes a renormalization-scheme parameter. 

The short-distance pressure pe is given 9 ) as an expansion in powers of the QCD 
couplant a s : 



PE 



PSB 



PE 



PSB 

8vr2 
45 

7T 
+ 



1 



15(1 + 5^/12) „ 
4(1 + 21^/32)^ 



21 

1+ 32 n/ 



rp 1 



(2-20) 
(2-21) 



1 + 



7T 



ttPq In 



+ 



7T 



7T 2 /3l 



(vr/?o) 2 ln 



72vr/3 



2vrT 

2 

2vrT 
1 + n f /6 



36 



1 + n//6 
l + 5n//12 

+ F(A E /(2irT)) 

Ae 



ln — + P h (n f ) 



1 + 5n//12 V 2vrT 



ln ^ + ^(n/) 



In 



2vrT 



+ S E 



(2-22) 



where 



F(A e /2ttT) 



244.9 + I7.24n f 



0.4150n2 



135(1 + n//6) 
478.4 + 35.14n/ +4.081n 2 f 



l + 5n//12 



In — + P h (n f ) 



(2-23) 



The renormalization-scheme parameters for pe is fi and /?2- Eq. (|2-22|> includes the 
unknown number 5e, which depends on the parameter p\. To see this dependence, 
we use the fact that the physical quantity pe, Eq. (|2-20|1 . is independent of the 
parameters, \i and /?2, up to and including 0(a^). Then from Eqs. ()2-5j) and (|2-22f) . 
we find the relation between 5e in the renormalization scheme with /?2 and its MS 
scheme counterpart (o~e)ms'- 

■2 



7T 

$E = (Se)ms --a-l&i- (A^MS 
A) 



(2-24) 



For the unknown (£e)ms> following Ref. 9), we will allow the following variation: 

-\F(Ae/(2ttT))\ < (5 E ) m < +\F{A e /{2t:T))\ . (2-25) 

While the pure number (Se)ms * s independent of any scale, F will have only slight 
dependence on temperature T. 

Finally, we note that, when the expression Pm + q, Eq. (|2-17|) . is used for the long- 
distance pressure, the total pressure Pm+g +PE is independent 9 ) of the "separation 
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(2) 

scale" Ae as it should be. It is worth remarking in passing that, if \i in p M+G is not 
equal to fi in pe, this is not the case. 

§3. Optimum approximants 

It is found in Ref. 6) that various approximants, which include the optimal 
approximant, constructed using the untruncated form (|2-llj) - (|2-12j) for pm+g are 
more reliable than those constructed using the truncated form (j2T7jl - (|2-19j) . It is 
found in Ref. 9) that this is also the case for the approximants constructed on the 
basis of the Pade and Borel-Pade resummations of pm+g- Keeping these observations 
and the remarks made above at the end of §2.1, we will also take more "reliable form" 
(HHU) - (j2~T2|) for pm+g, and apply the PMS criterion to it. 

Before entering into our analysis, it is worth summarizing here the procedure 
and the results of applying the PMS criterion in Ref. 6). 

1. The PMS criterion is applied to the total pressure p = pe + Pm+g- 

2. The 2-loop running couplant cr s (/i) is used, so that the renormalization-scheme 
parameter is [i only. 

3. Ae is chosen to be Ae = /i- 

The resultant optimal approximants are close to the lattice results for T > 0.7 GeV. 
The three-loop optimal pressure (with three quark flavors) is ~ 30% larger than the 
lattice result at T ~ 0.3 GeV, and the four-loop optimal pressure in pure-glue QCD 
with appropriate choice for Se and Sm+g is ~ 20% larger at T ~ 0.5 GeV. 

We will improve 1., 2., and 3. above as follows: 
1'. The PMS criterion is applied to the short-distance pressure pe and the long- 
distance pressure pm+g separately. 

(3) \ 

2'. For pe, the three-loop running couplant a s (/i) is used*- 1 , while, for pm+g, the 

(2) 

two- loop as {p) is used (see above after Eq. (|2-19l) ), 

3'. Ae is determined by requiring the minimal sensitivity of p on Ae (cf. Eq. ()3-4|) 
below). 

In the following,, the "long-distance" ("short-distance") /U and a s are denoted 
by hm+g and af +G (he and af), respectively. The OPT equations for pm+g read 
(see above after Eq. (|2-19j) ) 

dpu+G n daf +G (2) m+G\ 
HM+Gx =0) fJ-M+G^ =(3 (> {a s +u ), (3-1) 

O^M+G O^LM+G 



where Pm+g is as in Eq. (|2-llj) with Eq. I|2-12|) . The OPT equations for pe read 

dpE dp E . . 

VE^ — = 0, —— = 0, (3-2) 

dfj,E d(3 2 

*' As observed above in conjunction with Eq. 12-24H . when the PMS criterion is applied to pe 
(Eqs. (I2-2U> and 12-221 1 in a consistent manner, the 3-loop running couplant should be used. 
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where p E is as in Eq. (jlTlfflj) with Eqs. (jMlj) - J222J). 

The scale parameter Ae is arbitrarily chosen in the range 0{g s T) < O(Ae) < 
0(T), so that p = pm+g + Pe should not depend on Ae- As has been mentioned 

(2) 

at the end of §2, this is the case provided that (a) the truncated expansion p M+G 

1/2 

(in powers of a s ) is used for Pm+g an d He = Hm+g- I n our case, however, 
both (a) and (/3) are not fulfilled, and then p does depend on Ae- To resolve this 
ambiguity, we introduce a criterion for determining Ae- Respecting the spirit of the 
PMS criterion, we employ the following determining equation for Ae, 

Solving a set of equations Q3Tp - (j3-4j) . we obtain the optimized hm+g (= (J>l(T)) 



a 



at(ji L (T))), n E = Ms(n (3 2 = 02 S (T), af = of (fi S (T), /3 2S (T)) and 
A E {T) (= A L J S {T)). Here "L" and "S" stand for "long-distance" and "short- 
distance", respectively. Using thus determined quantities in Eqs. (|2T1|) and (|2-2Uj) . 
we obtain the optimized pressure p(T) = Pm+g(T) + pe(T). Since p = Pm+g + Pe 
is linear in Ae, when the hm+g , ^E and (3% are fixed to the respective optimized 
values, p(/j,l(T), /jis(T), @2s(T), Ae;T) is independent of Ae- 



§4. Results 



We are concerned with the temperature range 0.3 GeV < T < 3 GeV. Un- 
less otherwise stated, we will take for the number of active (massless) quark flavors 
nj = 3, which is of physical interest. For the QCD couplant a s (fi), we take the ref- 
erence value (a s (l.2 GeV))^jg = 0.403. The form for the short-distance af (he, fy) 
is obtained from Eq. (|2-7|l . which is give in Appendix A. Following the procedure 
stated in §3, we obtain two solutions within the reasonable ranges for hm+g, I^e 
and Ae- We refer the solution with larger (smaller) optimal approximant for p the 
solution I (II). 

For the time being, we choose 5g = (<5e)ms = Similar behaviors are observed 
for other values of 5q and (5e)m<t- 

In Fig. 1, we present the solutions I and II against T. The results are normalized 
by the Stefan-Boltzmann pressure psb (cf- Eq. ()2-21|) ). Also shown are the optimized 
long-distance pressure pm+g and the short-distance pressure pe- The circles are the 
predictions of the lattice calculations 11 ^ (see Ref. 9) for details.) From the figures we 
see that i) for the solution II, pm+g < and ii) the solution I gives better agreement 
with the lattice result. From these observations, in the following, we adopt the 
solution I as the better candidate for the optimal approximant. 

In Fig. 2 the optimized hm+g, Hl(T), and he, hs(T), are shown against T. 
The dot-dashed lines shows the traditional choice /u = 2-kT. The ratios /is(T)/T 
and hl{T)/T decrease monotonically with T: The ratio ns(T)/2irT is larger than 1 
at T = 0.3 GeV, ~ 1 at T = 0.5 GeV and becomes smaller than 1 for T > 0.5 GeV. 
Hl(T) is smaller than Hs(T) and he(T)/ ns(T) ~ 0.31 ~ 0.36 in the range shown in 
the figure. This is in accord with an intuition: The scale he that governs the physics 
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in the hard region, being of 0(T), is larger than the scale hg+m in the soft region. 

Fig. 3 depicts the optimized ft, 02s(T), against T. The dashed line shows 
@2 = ($2) ms- ftsC? 1 ) monotonically decreases with T crossing the MS value at 
T ~ 0.5 GeV. We note that Eq. (|2-5f) tells us that increase (decrease) in ft acts 
as increasing (decreasing) a s , the same effect as arising from decreasing (increasing) 
\i. Then, we suspect that, if we fix ft to (ft)Mg and vary only /z, we obtain much 
smaller (larger) optimized he for T < 0.5 (T > 0.5) GeV. Thus, we expect that the 
optimized he tends to 2irT. It is amazing to note that, at T ~ 0.5 GeV, Hs(T) and 
@2s(T) are close to the traditional choice \i = 2vrT and the MS value, respectively. 

L/S 

In Fig. 4, we show the optimized Ae, A J (T), and rriE{T) against T. For 

LIS 

reference Hl{T) in Fig. 2 is also shown again. We see that A E ' (T) 
Hl(T) (< hs(T))- This is not in accord with an intuitive expectation hm+g < < 
He- However, as seen at the end of §3, when hm+g, Ve, and ft are fixed to their 
respective optimized values, the pressure p does not depend on Ae- 

In Fig. 5, we show the pressure p at T = 0.5 GeV near the optimized values 
of of ~ a?(T), af + G ~ c^(T), ft ~ (3 2S (T) and A E * A L J S (T). In Fig. 5 (a), 
we show the behavior of p against a~ , when other quantities, a^ +G , P2 and Ae are 
fixed to the respective optimized values. The arrow indicates the optimized value of 
of. Similarly, Figs. 5 (b) and (c) show the dependencies of p on and on ft, 

respectively. From Fig. 5(c), we see that p depends on ft only weakly. As has been 
repeatedly pointed out, p does not depend on Ae, when af , a^ +G and ft are fixed 
to the respective optimized values. 

In Fig. 6, we present the optimal approximants for p, when the unknown con- 
stants (<5e)ms ana - are y aried. Within the ranges, Q2T6JI and (|2-25|) . of the 
unknown constants, (Sq, (^b)ms) = — F) yields the maximum p/psB, while 

(Sq, (<5e)ms) = ( — 5,+F) yields the minimum p/psb- For comparison, the MS re- 
sults are also shown*) . From the figure, we see as a bonus that the resultant optimal 
approximants depend rather weakly on Sg and (<5e)ms- ^ ur results are substantially 
improved when compared to the MS results. Taking into account that the lattice 
results should be taken within an error of 10 — 15 %, our results are well compatible 
with the lattice results. 

Finally, in Fig. 7, we depicts the optimal approximants (solutions I) for the case 
of pure glue QCD (rif = 0). For comparison, the solution II for 5q = (<5e)ms = i s 
also shown. Same comments as for the n/ = 3 case applies here. For 5g = (<5e0ms = 
0, the MS result accidentally coincides with the solution II with high precision. 
Although the results are substantially improved when compared to the MS results 
(except for the case of (So, (^-e)ms) = agreement with the lattice results 



pm+g (Eq. i'2-1'21 ) contains the logarithmic factors lnm_B and hig%, for which we have 
expanded as In a s + c\ + C2Ct s + ... with respective constants ci, C2 , ... . Then, for pm+g, we have 
used the truncated forms I2-17H - 12-191 1 and then p (= pm+g + Pe) is independent of Ae- For the 
f3{a s ) function, we have used the three-loop form. For the scale parameter /x, we have used the 
conventional one fi = 2nT. The unknown constant (<5-E)jyrg in Eq. 12-22L we vary in the range 
lEU. For A E in F(A e /2tyT) there, we have used Eq. | |2"3%1 with A E = ^/2-KTm E {T). For (j, M +a 
and a.s' I+G in ?ti_b(mm+g, a^ f+G ; T), we have substituted the respective optimized values. 
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are less than that for the realistic case of n/ = 3. 

§5. Discussion and Conclusion 

The QCD pressure p has been computed within the improved perturbation the- 
ory using the MS renormalization scheme, up to and including the 0{g® ln(l/g s )) 
and part of the 0(g®) terms. The truncated perturbation series for p does not rec- 
oncile with the relatively low T lattice results. As has been mentioned in §1, various 
approaches are proposed for improving the perturbative pressure p. Each approach 
performs a sort of resummation on the basis of respective philosophy. In this paper 
we add one more to them. 

We have applied the Principles of Minimal Sensitivity (PMS) criterion to p. 
This is not new, but in Ref. 6), among other things, the PMS criterion has already 
been applied to p. However, as has been mentioned in §3, application is not made 
in a consistent manner. Meanwhile, an important observation is made in Ref. 9): 
The contributions to p from the hard and soft energy regions, pe and Pm+g-, are 
separately physical quantities. On the basis of this fact, we have applied the PMS 
criterion to pe and pm+g separately. For determining the factorization scale Ae, 
which separates hard and soft scales, we also adopt the minimal sensitivity criterion. 

We have found that the optimal approximants are significantly improved when 
compared to the MS results. Especially, for the realistic case of rij = 3, agreement 
between our approximants and the lattice results is satisfactory. Moreover, we have 
found that, when compared to the MS results, the dependence of the approximants 
on the unknown terms of 0(g^.) is rather weak. On the other hand, for the case of 
nf = 0, although the results are substantially improved, agreement with the lattice 
results are poorer than that for the rif = 3 case. 

Further improvement may be expected if we perform the Pade type resummation 
starting from the optimal approximants obtained here, which is outside the scope of 
this paper. 
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Appendix A: Running Couplant 

The relation between the three-loop running couplants a s (/io, (/?2 ano - 

a s (/i, 02) is obtained by integrating the following two equations (cf. Eq. (|2-7|l ). 

da. 



/3 (3) («s), (A.l) 



W% = I2 («s) , (A.2) 

along a path C in a (ha.fi, p^-plane, which starts from the point (ln/xo, (/?2)ms) ano - 

arrives at the point (In /i, P2). Because of the integrability condition 80^ (a s )/d02 = 
(3) 

c?72 (a s ) /d In fi, one can use any path C. For convenience, we choose C = Ci©C2©C3 
withCi = (ha. no, (^)ms) (ln/u M , (P2)ms)> C % = ( ln MM, (A^ms) ~* ( ln A*M,/?2) and 
C3 = (ln/i a/, /%) — > (ln/i,/?2). At the final stage, we take the limit [iM +00. 
Integrating Eq. ()A.1|) along the path C\, we obtain 

fa l n Z' - 1 1 , A ,„ f «sM ^ /3l ^ f 00 + PlOLsM + (/?2)mS a sAf 



1__J_ + |i 1 /o £ m\_ A_ ln 
<sM a s0 Po V a s0 J 2 Po 



^0 7 a sM a s0 ' /3 *" V a s0 J 2/3 ™ V 0o + p\a s o + (P2)ms a s0 
;2 



-— ffl2^ MS [ G (^)ms) ~ G (°*0, (A)ms)] , 

z Po 

where a s0 = a s (jj, , (P2)ms) and = a * (W' (^)ms) and 



(A.3) 



2/?p / 2xy + 1 \ 

G(x,v) = — — arctan — — 

V^oy - M \y/*Poy-P'i) 

Integration of Eq. (jA. 1|) along the path C3 yields a similar equation for 0o ln(/i//ijw). 
Adding this to Eq. (|A.3|) . we obtain the formula for O ln(/z//z ), which includes, 
among others, a s ([j,m> (/^ms) an d oi s (iim,02)- From Eq. (|A.2I) . we see that 
a s (nM, (02)m) = a s{HM,02) + 0(a 3 s ). Substituting this for a s (fi M , (p^ms) in tne 
formula obtained above and taking the limit — ► +00 (a s (/j,M,02) — * 0), we see 
that the O(a^) term does not survive. 
Thus, we finally obtain 

a ln / M = J L + #L jjY _ 0\_ ln ^ A) + + AOs 



/io/ a s «so /So \asoy 2/3 \Po + /?i«so + (A2)ms a s0 



0j - 20 o (0 2 ) m 
01 ~ 20002 



[C(0, (02) m )-G(a sO ,(02) m )] 



2/3 2 

where a s = a s (/i, /? 2 ). 



[G(a s ,02)-G(O,0 2 )}, (A.4) 
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Fig. 1. The solutions I and II for the optimal approximant, p/pgg, as functions of T, when 
Sg = (<5s)ms = ^' (^SB i s the Stefan-Boltzmann pressure.) The optimized long-distance 
pressure Pm+g and the short-distance pressure ps are also separately shown. The circles are 
the predictions of the lattice calculations. 
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Fig. 2. The optimized ^m+g, Vl{T), and fi E , Hs{T), for S a = (fe)^jg 
dot-dashed lines shows the traditional choice fi = 27rT. 



as functions of T. The 
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Fig. 5. The pressure p near the optimized values of af, af I+G and /?2, when T = 0.5 Gev. The 
arrows in Figs, (a), (b) and (c) indicate the optimized values of af , af /+G and /?2, respectively. 
Details are given in the text. 
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Fig. 6. Optimal approximants (the solutions I) for p/psB as functions of T for various values 
of 5g and (fe)jyrg. Within the ranges, l|2-16[l and H2-25fl . of (8g, (Se)*^), (Sg, (fe)j^jg) = 
(+5, —F), yields the maximum p/pgg (dot-dashed line) while (Sg, — ( — 5, yields 

the minimum one (solid line). For comparison, the lattice results and the MS results are also 
shown. Displayed MS results are the maximum one [(Sg, (fc)fjg) — (+5, — -P 1 )]) within the 
ranges 12-161 and 12-251 . the minimum one [(Sg, (fe)}jg) = ( — 5, +F)] and the one with Sg = 




Fig. 7. Analogous to Fig. 6 but with nf =0. The solution II for Sg ~ (£e)tto = is also shown. 



